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We investigate the fast transport of an atom or a packet of atoms by different kinds of non¬ 
harmonic traps including power-law traps. The study is based on the reverse engineering method. 

Exact results are obtained and applied to design robust transport protocols. The optimization of 
the transport trajectory is performed with classical trajectories, and remains valid for the transport 
of a wave packet. 


I. INTRODUCTION 

The development of quantum information processing 
requires the accurate control of the motion of atoms 
or atomic packets [J-@. Similarly, cold atoms experi¬ 
ments often use the transport of atoms from a prepara¬ 
tion chamber to a science chamber [&, [tJ- To increase 
the number of quantum manipulations feasible in a given 
amount of time, a fast transport ending in a state without 
excitations is highly desirable. The experimental demon¬ 
strations of such diabatic transport have been realized 
with cold atoms transported in an optical tweezers [§} 
and for one and two ions with time-dependent electro¬ 
magnetic traps [ffl 

Recently, new proposed protocols to ensure the fast 
and optimal transport of neutral or charged particles 
have been workout. They are inferred from the differ¬ 
ent theoretical frameworks developed for shortcuts to 
adiabaticity [12] ■ A simple method is provided by the 
compensating-force approach fl2-fl6[. It requires to su¬ 
perimpose a time dependent constant force during the 
transport so to compensate exactly for the inertial force. 
The fast-forward formalism proposes the same solution 
E3- In practice, this trick may not be so easy to im¬ 
plement with the required accuracy. Alternatively, the 
reverse engineering approach, the Lewis-Riesenfeld in¬ 
variants scheme or the IV-point protocol provide a time- 
dependent tra p tr a jectory that ensures an optimal trans¬ 
port [H, [li], IT2 . fl9j . The direct implementation of 
the counteradiabatic approach [2Cj yields an unphysical 
Hamiltonian that can be recast as the reverse engineered 
Hamiltonian using an appropriate unitary transforma¬ 
tion [21.]. Those latter techniques have been developed 
for quantum and classical mechanics with a harmonic 
trap and for a one body problem. So far the role of 
anharmonicities has been investigated for trap expan¬ 
sion/compression in a Gaussian beam [22| and for the 
transport of two-Coulomb interacting particles [l5[ . 

In this article, we provide an extension of the reverse 
engineering approach when the transport is carried out 
with a non-harmonic trap. This study is important for 
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the experimental implementation since a fast transport 
implies that the high energy part of the potential is nec¬ 
essarily explored in the course of the transport while the 
harmonic approximation is only valid in the low energy 
range. We also investigate the robustness of the trans¬ 
port of a packet of atoms in the framework of classical 
and quantum mechanics. 

The paper is organized as follows. In Sec. [Hi we ex¬ 
tend the one-body reverse engineering method to non¬ 
harmonic traps, and develop an exact optimal scheme for 
a trap that has well-characterized anharmonicities. We 
then develop another protocol based on a perturbative 
approach about a fast and optimal trajectory for a har¬ 
monic confinement that is well adapted for a large range 
of anharmonicities. We also address the transport with 
a potential having a finite depth through the example 
of optical tweezers. In Sec. mu we study the robust¬ 
ness of the exact treatment developed in Sec. HH for the 
transport of a packet of atoms both with a harmonic and 
non-harmonic trap. In the last section, we show that 
the previous schemes based on classical mechanics are 
still valid for the transport of a wave packet in a non¬ 
harmonic trap. 

II. ONE-BODY TRANSPORT IN CLASSICAL 
PHYSICS 

A. Reverse engineering protocol with power-law 
traps 

Consider the transport of a particle of mass m confined 
in a power-law trap with an even exponent: U n (x) = 
mr] n [x — xoit)] 2n /{2n) with n > 1 an integer where Xq (t) 
denotes the trajectory of the bottom of the trap to be 
determined. We consider the transport over a distance 
d in a time interval tf. According to Newton’s law, the 
particle obeys the differential equation: 

x + r/n [x- xoit)} 2 ' * 1 - 1 = 0 . ( 1 ) 

The boundary conditions for the bottom trap trajectory 
are £o(0) = 0 and xo(tf) = d. To ensure an optimal 
transport, the particle should be at rest at t = 0 and 
t = tf, and obeys therefore the boundary conditions 
x(0) = x(0) = x(0) = 0, x(0) = d, x(tf) = x(tf) = 0. 
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FIG. 1. Robustness of the transport in quartic potential (n = 
2) against fluctuations of the final time. The relative error 
\x(t)/x(tf) — l\ is plotted as a function of t/tf for a polynomial 
interpolation of order 5 (solid line), 7 (dashed line) and 9 
(dotted line). 


The reverse engineering protocol works as follows. We 
set the trajectory x(t ) of the particle according to the 
boundary conditions using a time interpolation. For sake 
of simplicity, we take a polynomial function: 

x{t) = d [m/t f f - 15 (t/tf)* + 6 (t/tf) 5 ] . (2) 

We then infer the bottom trap trajectory: 

fex\ 1/(2n - 1] 

x 0 (t)=x(t)+e — , (3) 

VW 

with e = +1 for 0 < t < tf/2 and e = —1 for tf/2<t< 

tf- 

The robustness against final time fluctuations can be 
improved with a polynomial interpolation of higher or¬ 
der associated with the cancellation of the next order 
of the derivative of x at t = 0 and t = tf. For in¬ 
stance, for a polynomial interpolation of order 7, the 
extra boundary conditions, 'i'(O) = x[tf) = 0 are in¬ 
cluded. This is illustrated in Fig. [T] for the quartic po¬ 
tential n = 2. The expansion about the final time tf 
scales as \x(t)/d — 1| ~ \t/tf — l\ p ~ 2 where pis the degree 
of the polynomial that is used for the interpolation. How¬ 
ever, according to our numerical simulation this method 
does not improve the robustness against fluctuations of 
the trap strength rj n . 


B. Anharmonicities 

Consider now another type of non-harmonic potential 
that is the sum of a harmonic confinement and a cubic 
anharmonicity: 

U (x) = ^rnujg (x - x 0 (t)f + (x - x 0 (t)) 3 , (4) 


where £ quantifies the strength of the anharmonicity. 
Such a potential can result from an expansion around 
the minimum of the real transport potential (23] . The 
equation of motion reads 

w 2 

x +U}Q (x - X 0 {t)) + -J- (x - x 0 (t )) 2 = 0. (5) 

From this expression, we can infer the value of xq (t) if 
we use the polynomial interpolation ©) for x(t). This 
equation is of the form X 2 + £X + q = 0 with X = x — Xq 
and q = £ x/lUq. The solution reads: 


xo (t) 



( 6 ) 


This solution exists only when the discriminate is posi¬ 
tive i.e. for t;ujQt 2 f/d > dO/v^ according to the polyno¬ 
mial ansatz 0- This inequality reflects the fact that the 
harmonic plus cubic potential has a finite depth. The 
acceleration that can be used for transport is therefore 
bounded. In the limit £ —> oo, Eq. © collapses to Eq. © 
as expected. 

Consider now that the anharmonicity is quartic. We 
need to add to the harmonic potential a potential of the 
form mwg (x — xq( t)) /(4£ 2 ). Repeating the previous ar¬ 
gument, we obtain the bottom trap trajectory by solving 
a cubic equation: 


x 0 (t) = x(t) + 


£2/3 


2 1 /3 w 2/3 


f 2 , ,4 i 07 "2 \ 1/2 j ^ 


4£ 2 wg + 27ar 

W 


£2/3 


2 V 3 o; 2/3 


4£ 2 Wq + 27x 2 
27~ 


1/3 


■ (7) 


As previously, in the limit £ —> oo , we recover Eq. ©• 
The solutions presented above assumed that one knows 
exactly the value of the anharmonicity. In this case, we 
have shown that an exact strategy can be worked out us¬ 
ing the reverse engineering approach. However, one may 
need an approach that works at best but only approxi¬ 
mately for a large range of anharmonicities. We explain 
hereafter how one can proceed in this latter case. 


C. Perturbative approach to overcome 
anharmonicities 

To minimize the sensitivity to the anharmonicities for 
the final position and velocity of the particle, we use a 
perturbative expansion about a fast and optimal strategy 
designed for a harmonic confinement. We work hereafter 
out this method explicitly for cubic anhamonicities. The 
results of a similar treatment applied to the case of quar¬ 
tic anharmonicities are then given. 

First, we use as a reference the 1-point protocol of 
Ref. pjl to design the trajectory for transport in a har¬ 
monic trap of angular frequency wq. In this method, one 
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defines an auxiliary function that obeys the boundary 
conditions g( 0) = g(tf) = g{ 0) = g{tf) = 0 and the 
integral relations 

[ g(t)dt = 0, and [ At' [ g(t")dt" = -i. (8) 

Jo Jo Jo uo 


The 1-point protocol gives the trajectory through the 
differential equation xq (t) = g + to^g. By construction, 
xo(t) obeys the boundary conditions: xo(0) = 0 and 
xo(tf) = d. Using the simple polynomial interpolation 
g(t) = — t/tf) 2 { 1 — 2 t/tf) with J\f = d/ujQ A 

and A = 1/420, we obtain the following polynomial form 
for the trajectory of the bottom of the trap: 


x 0 (s) 


420 



-4s 3 + 





with s = t/tf, Xq(s) = Xo(t)/d and u = u>otf. With this 
choice, the trajectory of the particle can be explicitly 
worked out X\ (s) = 35s 4 — 84s 5 + 70s 6 — 20s 7 . One 
can check that X\ + Wq(xi — Xo(s)) = 0 with Xi(0) = 
0, Xi(0) = 0, Xi(0) = 0, xi(l) = 1, Xi(l) = 0 and 
xi(l) = 0 (the double dot corresponds here to a second 
derivative with respect to the s variable). We then solve 
perturbatively Eq. ©: 


~ 2 /~ _ . uj^d,_ _ . 2 d (xi) 2 

x 2 + u 2 (x 2 - x 0 ) =-Mx 2 - x 0 ) - -7— 2~- 

£ £ u 2 

The perturbative solution to the first order is x 2 (s) = 
$i(s) + (d/£)/i( s ) with 

/i(s) = —^ [ x 2 (s') sin[u(s — s')]ds'. 
u Jo 

The quality of the transport is evaluated by calculating 
the residual energy communicated by the transport at 
the final time f/, 


A E mtuod 2 x 2 

tvjjQ h 


(x — x 0 ) 2 d 


The optimal strategy that we propose proceeds in the 
following manner. First, we compute the dimensionless 
quantity AE/Hloq for the perturbative solution x 2 for dif¬ 
ferent values of £. Second, we search for the optimal val¬ 
ues, u 0 , of u = xiotf that minimizes the excess A E of 
energy after the transport. As illustrated in Fig. [2] (a), 
in the time interval 0 < tf < 4 x (27r/wo), we find two 
values of the final time tf ~ 6.97/wo and tf ~ 21.21/wo 
for which the transport is ultra robust against the an- 
harmonicity parameter £. They corresponds to a highly 
non-adiabatic transport. For the optimal choices uq the 
final excess of energy scales as A E ~ (d/t;) 4 hu! 0 since the 
lowest order contribution of the anharmonicities ( d/i.;) 2 
is cancelled out by the choice u = uq. We conclude that 
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FIG. 2. Residual energy, A E, after transport normalized to 
hcoo as a function of the anharmonicity parameter £/d and 
for different transport time tf. (a) cubic anharmonicity and 
(b) quartic anharmonicity. The calculations are performed in 
dimensionless units and are therefore generically valid. The 
prefactor for the quantitative estimate of the excess of energy 
have been estimated with the following parameters loo = 2nx 
1.41x10 s , m = 40xl.667xl0“ 27 ( 40 Ca+), a 0 = [h/(muJo)} 1/2 
and d = 20.2 x ao- 


the optimal non-linear strategy that we have derived pro¬ 
vides a set of discrete optimal transport times. For a 
longer transport time compared to the optimal one, ex¬ 
cited states can be massively populated. Fig. [2] (a) also 
provides a window of final time values about the opti¬ 
mal one for which the transport remains quite optimal. 
The low values of £/d correspond to large anharmonici¬ 
ties. However for a cubic anharmonicity, one has only a 
limited range of possible values since the potential expe¬ 
rienced by the atoms should remain a trapping potential 
in the course of the transport. 

From this respect, the situation is quite different with 
a quartic anharmonicity. The robustness of our approach 
can be tested outside the perturbative regime. Repeating 
the previous argument for the quartic anharmonicity, the 
lowest order perturbative solution reads x 2 (s) = Xi(s) + 
( d/ 0 2 f 2 {s ) with / 2 (s) = (1/w 5 )/p (xi(s')) 3 sui[m(s — 
s'jjds'. The optimal value uq that minimizes the excess 
of energy at the final time 


A E 
hu> 0 


mujod 2 

h 


x 2 

~2u 2 
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after the transport is uq — 14 in the time interval 
0 < tf < 4 x (27 t/wo) (see Fig. [2] (b)). The excess of en¬ 
ergy, AE/Hujq, decays as (d/C) 8 about the optimal value 
while the lowest order contribution of the anharmonic- 
ities provides a scaling as (d/C) 4 - Figure [2] (b) shows 
explicitly the breakdown of our approach for too large 
anharmonicity (C < 0.3d). Interestingly, new features 
outside the perturbative regime emerge as nearly hori¬ 
zontal lines. They are associated with specific values of 
the anharmonicity (for instance C = 10~ 173 d) for which 
our designed trajectory provides an optimal transport ex¬ 
tremely robust against the final time (A E < 10 _11 /iWo 
for 15.5 < ujotf < 21.6). 

The method developed here can in principle be further 
improved by searching for a solution of the next order in 
the expansion of x 2 (s). The optimization procedure re¬ 
quires to choose the final time u and the other parameters 
to cancel higher order terms in the expression for the ex¬ 
cess of energy after the transport. For this purpose, one 
needs to add extra free parameters in the g function that 
generates the xq trajectory. Alternatively, one could ap¬ 
ply the strategy for optimizing the spring-constant error 
in Ref. f2dl , in which the polynomial form of higher order 
is designed to nullify the integral /i or / 2 . As a result, 
the residual energy A E can be reduced by one order of 
the magnitude, as compared to the original polynomial 
form, Eq. ©■ Finally, the ultra robustness against final 
time observed outside the perturbative regime for quar- 
tic potential (Fig. [2] (b)) can be tuned by adding extra 
parameters to the trajectory on which the protocol is 
based. 


As previously, we shall use the polynomial interpolation 
© for x(t). To extract the optimal trajectory using the 
reverse engineering method, we first consider the time 
interval 0 < t < tf /2 for which x{t) > 0. One can 
readily show that Eq. m has solution only if {d/t 2 ) < 
(27-\/3/1040)77. Physically, this latter criterium sets an 
upper bound for the ratio d/t 2 resulting from the finite 
depth of the potential. Indeed, this potential remains a 
transport potential only for acceleration below 2777/104 
according to our ansatz ©. In the appropriate range 
of acceleration, the trajectory is obtained by solving the 
quartic equation m in X. The analytic solution can be 
readily worked out, and the time evolution of the focal 
point is then given by Xq (t) = x(t) — X(t). 


III. TRANSPORT OF A PACKET 

In this section, we evaluate the quality of the fast trans¬ 
port protocol for a packet of atoms using the reverse engi¬ 
neered protocol designed in the previous sections. In the 
first subsection, we explain the new features that arise 
when one consider the transport of a packet of atoms 
using a harmonic confining potential, we then present 
a strategy to optimize the transport in the presence of 
anharmonicities. 


A. The transport of a packet in a harmonic trap 


D. Transport in optical tweezers 


The first transport experiment of cold atoms carried 
out outside the adiabatic regime was performed with 
an optical tweezers generated by a focus Gaussian laser 
beam whose focal point was displaced using an accurate 
translation stage [8f. The potential experienced by the 
atoms along the longitudinal axis is 


U(x) 



Up 

[x - X 0 (t )} 2 ’ 


( 10 ) 


where xr = nw^/X is the Rayleigh length of the Gaus¬ 
sian beam, wo is waist and A its wavelength. This poten¬ 
tial is clearly non-harmonic and has also a finite depth. 
The three first terms of the expansion of the potential 
m around its minimum provide a quartic potential for 
which the previous analysis can be used. Alternatively, 
the transport can be solved exactly in this specific case. 
Introducing the variable X = [x — Xq( t)\/xR and the pa¬ 
rameter rj = ZUo/mxR, the equation of motion reads: 




( 11 ) 


To evaluate the quality of the transport we shall com¬ 
pare before and after the transport the variance, (Ax) 2 = 
x 2 — x 2 , with 

1 N i N 

X 2 = X i and X = jy Xi ^>- ( 12 ) 

v i=l V i= 1 

First, we calculate exactly those moments as a function 
of time in the case of a harmonic trap of angular fre¬ 
quency lo o- The initial standard deviations of position, 
Axi, and velocity, Axj, are related at equilibrium. For 
instance, in a harmonic trap at thermal equilibrium, one 
has A Vi = UqA Xi = (fc^T/m) 1 / 2 where T is the temper¬ 
ature and m the atomic mass. An optimal transport for a 
packet needs two requirements: the center of mass should 
follow the optimal one-particle trajectory (see Sec. m, 
and the standard deviation should obey a relationship 
that characterizes the equilibrium after the transport. In 
the thermal case that we consider, we expect the equality 
Avf = wo Ax/. 

In the following, we propose a moment method that 
shows explicitly how the center of mass motion (x and v 
moments) is coupled to the quadratic moments (x 2 , v 2 
and lev) in the course of the transport. Starting with the 
calculation of the time derivative of the moment x 2 , we 
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derive the following close set of equations 


x 2 = 2xv , (13) 

XV = V 2 — U) qX 2 + OJqXoX, (14) 

v 2 = — 2cjqXV + 2loqX 0 v , (15) 

x = v, (16) 

v =—luqX + lvqX q. (17) 


Imposing that the packet is at equilibrium at t = 0 and 
at t = tf amounts to setting the values of the position- 
velocity correlation moment, xv, at the boundaries of 
the time interval (xU(0) = 0, xv(tf) = 0, xU(0) = 0 and 
xv (t f ) = 0) in addition to the boundary conditions on the 
bottom of the trap variable Xq( t). As a result, we recover 
from Eq. (fl4l) the equilibrium condition that relates the 
standard deviation for position and velocity at t = 0 and 
t = tf. v 2 (0) = UqX 2 (Q) and v 2 (tf) = ujQX 2 (tf). 

From Eqs. cam we deduce the equation for the 
position-velocity correlation moment: 


XV + AuJqXV = Wg (xqx + 2 >xqv) 


(18) 


whose exact solution is 

nt 

xv(t) = ^ / (xox + 3xov) sm[2u>o(t — t')]dt', (19) 

2 Jo 

where we have taken into account the boundary condi¬ 
tions SU(0) = 0 and auT(O) = 0. Interestingly, the poly¬ 
nomial ansatz m for x(t) provides a bottom trajectory 
xo(t) = x(t) + x/lOq such that the extra boundary con¬ 
ditions xv (tf) = 0 and xv(tf) = 0 are automatically 
fulfilled. 

The solution found here using a classical formalism is 
actually also valid quantum mechanically for the trans¬ 
port of any eigenstate of the harmonic oscillator. Indeed, 
the set of Eqs. (I13H17D can be derived using the Ehren- 
fest theorem for the observables X , P 2 , XP, X and P 
f25j , and the relation between the quadratic position and 
velocity v 2 = cvqX 2 remains valid for all eigenstates. The 
fact that an optimal choice for the bottom of the trap 
ensures the perfect transport of the position of velocity 
dispersions is a result specific to harmonic trapping and 
that can be proved using alternative approaches such as 
the scaling method (26j. This latter technique enables 
one to show that the design of the optimal transport for 
a one-body wave function is still valid for an interacting 
Bose-Einstein condensate in the Thomas-Fermi regime 

sa 


B. Transport of a packet in the presence of 
anharmonicities 

In this subsection, we use the exact results of Sec. Ill Bl 
for cubic and quartic anharmonicities for the bottom trap 
trajectory and study the quality of those optimal strat¬ 
egy for the transport of a packet of atoms. Our bottom 



FIG. 3. Relative excess of final energy after the transport as 
a function of the transport duration tf and the strength of 
the anharmonicity (a) cubic anharmonicity and (b) quar¬ 
tic anharmonicity. Dark vertical lines correspond to “magic 
times” for which the transport of the packet is optimal. 


trajectory choice provides a perfect transport of the cen¬ 
ter of mass of the packet. To evaluate the impact of this 
transport on a packet we calculate the relative variation 
|1 — Ef/Ei\ of the total energy after the transport com¬ 
pared to its initial value. This quantity is studied as a 
function of the strength of the anharmonicity and the 
transport time. The transport time plays the role here 
of a free parameter that will be adjusted to guarantee an 
optimal transport. 

The results are summarized on Fig. [3] (a) for the cu¬ 
bic anharmonicity and Fig. [3] (b) for the quartic one. In 
both cases, we observe the emergence of a discrete set of 
transport times that ensure the perfect transport of the 
packet with a very small excess of energy after transport. 
For example at £/d = 10 1 ’ 5 in Fig. [3] (a) , we find the 
following set of magic times («o,Log 10 |l — Ef/Ei\)) = 
(5.7635,-5.1974), (7.8343,-5.4588), (9.0968,-5.666), 
(10.995, —5.9142), (12.329,-6.0903), (14.1238,-7.2799), 
(15.513, —6.4508), (17.285,-6.8358), (18.69,-6.7855), 
(20.412,-7.2077), (21.8453,-7.793), (23.575, —7.1666). 
Those “magic times” are quite robust against the 
strength of the anharmonicity (this is the reason why 
they appear as vertical dark lines in Fig. 0 . The general 
conclusions deduced here in the specific case of cubic an¬ 
harmonicities are generic and holds in the case of quartic 
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FIG. 4. Transport of a quantum wave packet. Time evolution 
of the total energy, E(t), with respect to the initial energy, Ei, 
in the course of the transport for three different time transport 
tf = 0.75 t*f (dotted line), tf =t*f = 13.335/u;o (solid line) and 
t = 1.25t^ (dashed line) where tf has been obtained from the 
results of Sec. Ill Cl based on classical mechanics. 


anharmonicities as illustrated in Fig. [3] (b). 


IV. TRANSPORT OF A PACKET: QUANTUM 
ANALYSIS 

So far, the analysis has been carried out using classical 
mechanics. The global picture of the optimal trajectories 
depending on the different parameters could be obtained 
rapidly in this manner. An important question is the 
validity of this classical approach to transport a quan¬ 
tum wave packet. It would be too time consuming to 
run the 500 000 numerical simulations that have been 
necessary to realize the 2D plots of Figs. [2] and [3] We 
have checked on a few examples that the transport of a 
quantum packet is in perfect agreement with the classical 
prediction. An example is provided in Fig. [4] for which 
a relatively large quartic anharmonicity has been used 


(Logio(£/d) = —0.8). We have plotted the total energy, 
E(t), (relative to its initial value) as a function of time in 
the course of the transport for three different time trans¬ 
port tf = 0.75 t*j, tf = tf and t = 1.25 t*f where tj co¬ 
incides with the prediction of magic times deduced from 
classical mechanics (see previous section). The number 
of excited state transiently populated rangers from 4 to 
10 depending on the transport time and clearly shows 
the non-adiabatic character of the transport. The opti¬ 
mal value that we find for the transport time is neither 
the larger one nor the one that transiently populates the 
minimum number of excited states. This is a general 
feature of shortcuts to adiabaticity results. 

V. CONCLUSION 

In this article, we have investigated the use of the re¬ 
verse engineering method to perform the fast and ro¬ 
bust transport of an atom or a packet of atoms in the 
presence of non anharmonicities or with power-law trap. 
The results provide a clear strategy to optimize the non- 
adiabatic transport of atoms. Interestingly, the approach 
based on classical mechanics provides also the optimal 
solution in the quantum case for the fast transport of a 
wave packet. This work should be useful for the commu¬ 
nity of cold atoms manipulation and also for quantum 
information operations in which processing and storing 
sites are separated as in the case of ion trapping physics. 
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